Packing points into a unit cube in higher space. 
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Abstract 

In this paper using the concept of the extended Hamming code we give a construction 
for dense packing of points at distance at least one in such unit cubes which dimension are 
a power of two. 
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1 Introduction 

> 

. The fohowing problem was stated in [13] and later repeated in [S], |14j . [1]. 

■ Let f{n) denote the maximum number of points that can he arranged in the n- dimensional 

i unit cube (n-cube) so that all mutual distances are at least 1. Obviously, f{n) = 2" for n < 3. 

Many have shown that logf{n) ~ ^n{logn). Determine the exact values of f{n) at least for 

• small n. 

O ■ 

Any construction of a suitable point-set gives a lower bound for f{n). Previously constructed 
_ sets (see [3], [5] and [2], p.71) shows that /(4) = 17, /(5) > 34, /(6) > 76, /(7) > 184, 

^ '. /(8) > 481, /(9) > 994, /(lO) > 2452 /(II) > 5464, /(12) > 14705. To show a good upper 

d ' bound is usually much more difficult. In the paper [6] the authors proved the following upper 

estimates: /(6) < 192, /(7) < 576, /(8) < 2592, /(9) < 11664, /(lO) < 46656, /(II) < 248832, 
/(12) < 944784. More results are known for the dimension 5, for example /(5) < 44 (see ^). 
Better upper estimate /(5) < 43 was shown a short time ago in [9] and this was most recently 
strengthened to /(5) < 42 in [10]. The most recent result (see [7j) says that /(5) < 40. 

The best known asymptotic estimates can be found in [5] they are /(n) < n^/^O, 63901"e°(") 
and f{n) > n"/^0, 2419707"O(-y/n). The lower bound is not a constructive one. 

In this paper using the extended Hamming codes we give a construction for packing 

3" + 2(n - l)3"/2 + 1 
2n 

points at distance at least one in the n = 2'^-dimensional unit cube. In the cases of A; = 2, 3 the 
corresponding packing give the known better ones. This construction can be generalized by the 
Reed-Muller (briefly R-M) codes. Since the weight-distributions does not known in general, we 
give lower bound for the number of points of this construction. For odd k it is equal to 

3" + 2(n - 1)3-/2 + 1 ^ ^ i(i,g^„+i) 2 



+ / • n2^'°S2«+ii _ n[{n - l)(n - 2) - 3], 



2n 3 
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where I = 4.768462..., and for even k is 

3" + 2(n-l)3"/^ + l ^ _ 2) . „i(iog,n+i) _ 2 _ ^^^^ _ 2) - 3] + i^^. 

2 Constructions using the binary extended Hamming code 

The codes that Hamming devised, the single-error-correcting binary Hamming codes and their 
single-error-correcting, double-error-detecting extended versions marked the beginning of coding 
theory. These codes remain important to this day, for theoretical and practical reasons as well 
as historical. In our paper we give a new application of the extended version of this code, we use 
it to define large point system in the unit cube at distance at least one. To the calculation of the 
number of points wc have to use the weight distribution of this code, which is the distribution 
of the Hamming weights of the codewords belong to this code. So we recall to its generator 
function. 

A simple definition of the binary Hamming code is the following one. Let k be a positive 
integer and construct a binary k x (2^ — 1) matrix H such that each nonzero binary /c-tuplc occurs 
exactly once as a column of H. The codewords (of lengths 2^^ — 1) of the Hamming code are the 
binary — l)-tuples as column vectors orthogonal to the rows of H. It is easy to see that binary 
Hamming code is a linear code having minimum Hamming weight and distance 3, meaning that 
every codewords contains at least three nonzero coordinates. If we add an overall parity check 
bit to a binary Hamming code then we have an extended Hamming code. This means that for 
a codeword we glue a 2'^-th coordinate it is zero if the weight of the codeword is even and one 
if it is odd, respectively. The extended binary Hamming code is a mod2 linear vertex set of the 
n-dimensional unit cube each of which contains at least four nonzero coordinates. 

First we consider those points of the cube which are on the middle of the position vectors 
occurs as a codeword of the extended Hamming code. If the Hamming weight of a codeword 
is j then the corresponding point there is j coordinates with value \ and its other coordinates 
are zero. Corresponding to this codeword wc collect those points which other coordinates are 
either zero or one. Such a way a codeword with weight j generates 2 ~^ points in the unit cube. 
The distances of these points to each other at least one, since the minimum weight of the code 
is 4 showing that the substraction of two position vectors contains at least 4 coordinates with 
absolute value \ or at least one with absolute value 1. We now determine the number of points 
of this system. 

The weight for the Hamming code are found by use of the correspondence between linear 
dependence relations among r columns of the matrix H and codewords of weight r. The number 
of codewords of weight j is denoted Wij). The number of codewords that are zero linear 
combinations of (j — 2) vectors is W{j — 2). The number of linear combinations consisting of 
j — 2 terms that add to plus one nonzero term (n — (j — 2))W{j — 2) where n = 2'^ — 1. 
The number of linear combinations of {j — 1) vectors is The number of nonzero linear 

combinations of {j — 1) vectors is {j^i) — W{j — 1). The number of nonzero linear combinations 
of {j — 1) vectors to which one more term may be added to form a linear dependence among j 
vectors is therefore 

jW{j) = ^ J - WU - 1) - (n - (i - 2))W{j - 2). 

This recurrence relation gives a possibility to calculate the weights of this code. Of course, 
W{0) = 1, W{1) = W{2) = 0. To an explicit solution we introduce the generator function 

n 

f{x) = J2w{j)x^. 
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As it can be seen in [12j f(x) can be determined and it is 



f(x) = (1 + x)" + n(l + a;)("'i)/2(l - x)("+^)/2 

n + 1 L 

where n = 2^ — 1. 

Now the weights for the extended binary codes can be found by nothing that each codewords 
of odd weight has a 1 added to it, while each of even weight has a added. Thus the number of 
codewords of weight j is if j is odd and V{j) := W{j)+W{j — 1) if j is even. The odd terms and 
even terms can be separated by using [f{x) + f{—x)]/2 for the even terms and [f{x) — f{—x)]/2 
for the odd terms. Then the generator function as we can see in [15j is 

9ix) = lifix) + f{-x)] + |[/(x) - f{-x)] = 



1 



:i + x)''^ + (1 - xf + 2(2^= - 1)(1 - x^^^'~ 



2A:H 

The number of our points can be calculated on the base of the numbers V{j) since two points 
are distinct either they associated to distinct codewords or they are at least one distinct integer 
coordinate. Thus if we denote by H{j) the number of points associated to codewords with 
weight j we have: 

H{j) = 2^'-^Y{3). 



Equivalently, we have 



2«: 2^ 



j=0 j=0 

2"' g{x) = Y,H{mxy . 

j=0 



2k (l\ 3" + 2(n - l)3"/2 + 1 



2 J 2n 



Thus 

Y.H{j) = 2^^g 
j=0 

if n denotes the dimension 2^ of our space. 

This construction usable if k > 2 and gives point systems with cardinality 17, 481, 1351361 
in the respective dimensions 4,8 and 16. The first two values give the best known results, 
respectively, but in the third case we can give better point system. Observe that in a 16- 
dimensional space certain points which all coordinates are either ^ or | can be attached to the 
previous point system, such that, the distances of it from the original points are at least one. 
To give a relatively large system with at least one pairwise distances, consider the extended 
Hamming code which is a linear code with minimal weight four, and adding the the points of 
forms jv + j{l — v) to the system where u is a codeword of this extended Haming code. Since 
the subtraction of two codewords contains at least four nonzero coordinates the subtraction of 
such two point there are at least four coordinates with absolute value greater or equal to |, 
showing that the distance between the points is at least one. Thus it can be attached 

5(1) = 5^ y(i) = 22^-4-1 = 2048 
j=o 

new points and we get a better point system with cardinality 1353409. 
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3 A generalization of the construction 



As we observe in the previous paragraph, to the codewords of greater weights we can correspond 
to points containing smaher coordinates. So respectively handing the codewords with the same 
weights is a natural possibility. However the difference of two codewords with the same weights 
can be codeword with smaller weight, it is more clearly to use linear sub-codes of the Hamming 
code as the set of codewords with fixed weights. Fortunately in the extended Hamming code 
EH{2^) we can find nested sequences of linear codes with increasing minimal distances, 4, 16, . . . 
these are the Reed-MuUer codes denoted by RM[2'', 1 + Y.l=i for r = (A; - 2), (A; - 

4)...,fc-2[|]. 

In our construction firstly we collect the vertices of the original cube as points. Then we 
determine the further points of the arrangement corresponding to the element of the R-M codes 
mentioned above. Let v G RM[2^, 1 + YlUi it) ' 2*'~1 for an r G {{k - 2), {k - 4), . . . ,k - [|]}, 
the weight of v is w{v) / 0. To this codeword we correspond to all points of form 

+ (1 - v) ■ e, 



2^ 

where 1 is the codeword (!,...,!) of length 2*^, {l — v)-e means an element of the universe code 
of length 2^ — w{v) positioning to the nonzero coordinates of the codeword {1 — v), I is odd and 
1 < I < 2^^*"^^. Remark that a binary vector v will be take into consideration with respect to 
all of the codes containing it. The number of such points is 2~ ^ • 2^'°~"'(''), and in every steps 
we are corresponding to v such points in the unit cube which were not taking before. We also 
remark that the present construction is a simplified variation of that one which was investigated 
in the preceding paragraph. The reason is that the definition of the point system in such a way 
is more clearly. Now the total number of points in this construction is 



3" + 2(n - l)3'*/2 + 1 ^ ^ 



2n 



+ 



E 2^ 



E 

Vr=fc-2[|]/ "'"^^V^x 

\ ^ veRMl2k,l+Y:V^^(^),2k-r 



^-1 . r)2''-wiv) 



where means that we take the sum only for even or only odd integers, respectively. Un- 

fortunately there is no general formula for the weight distribution of the Reed-Muller codes 
except in the cases of r = 0, 1, 2 and k — 2. These are the repetition, even weight, second order 
Reed-Muller and Hamming codes, respectively. Thus there is no chance to determine the exact 
value of this sum. An asymptotic lower bound we can get using the minimal weight elements 
only since their number are known. Let A2k-r be the number of codewords of minimal weight. 
Then we have: 



A2k—r — 2 



k-r-1 r,k- 



2«-* _ I 



2k—r—i J ' 

i=0 

and thus a lower bound of the second part of our sum is 



Introducing the notation r' = k — r we get that this sum is: 



2'" J- -L 2^ 2* 

r'=4/ \ i=0 



For odd k the maximal value of r' is 2[|] = k — 1 meaning that if k > r' > 2 then 



L2J 

/ f? I 

2*^-2^ + !>_ + _ + 1 

2 ~ 2 2 



and thus the number of points is greater than 

3" + 2(n - l)3"/2 + 1 ^ ^4+^+1 / 2^ 'Jj' 2^= - 2 



2n \ ^ / \ 2*^ J-J- 2'- - 2* / 

\r'=4/ \ i=0 / 

3" + 2(n - l)3"/2 + 1 /^\* /2 /2'= '^-TT^ 2^ -2*\ 2 , , 

= + [gj (^2^ n j - - - 2) - 3] = 

= 3" + 2(n- 1)3-/^ + 1 ^ , , „„.^„,„ _ 2 2) _ 3,, 

2n 3 
where / = 4.768462..., as John Leech enumerated in 

Analogously it can be seen, that the lower bound in the case when the number k is even a 
little bit worth, it is 

3n + 2(n-l)3"/2 + i _ ^i(iog,„+i) _ 2 _ ^^(^ _ 2) _ 3] _ 2(^|(iog2 n+D _ i^) 

showing that the number of points for large n is greater then 

3" + 2(n - l)3"/2 + 1 ifw r,-4-i^ 2 , , I ^ 

^ — + - 2) • n2('°S2n+i) _ ^[(^ _ _ 2) - 3] + 

as we stated. 
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